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Abstract 

We introduce a new minimisation principle for Poisson equation us¬ 
ing two variables: the solution and the gradient of the solution. This 
principle allows us to use any conforming finite element spaces for both 
variables, where the finite element spaces do not need to satisfy a so-called 
inf-sup condition. A numerical example demonstrates the superiority of 
the approach. 
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1 Introduction 


It is often more important to get the accurate approximation of the gradient 
of the solution of a Poisson equation. In that case, a mixed formulation of the 
Poisson equation is used, where there are two unknowns - the solution and its 
gradient - in the variational equation. Discretising a mixed formulation of a 
partial differential equation is a challenging task as the involved finite element 
spaces should satisfy a compatibility condition - so called inf-sup condition. 
Although there are many finite element spaces discovered satisfying the com¬ 
patibility condition for the Poisson equation [2ll6lM [1^116] . it is not so easy 
for mixed formulations of other partial differential equations. It is sometimes 
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useful to use a least-squares finite element method to approximate the solu¬ 
tion and its gradient simultaneously [SHH]. A least-squares formulation allows 
the use of any conforming finite element spaces avoiding the compatibility 
condition. 

In this paper, we propose a new minimisation principle for Poisson equation 
using the solution and the gradient of the solution as two unknowns. This 
formulation is similar to a least-squares formulation in the sense that it allows 
the use of any conforming finite element spaces avoiding the compatibility 
condition [HHSlin]. However, in comparison to a least-squares finite element 
method, the source term / can be in the dual of i^^-space, and the gradient 
can be discretised using a L^-conforming finite element space. We also prove 
optimal a priori error estimates for the proposed finite element method. 

The structure of the rest of the paper is organised as follows. In the next 
section, we introduce our formulation and show its well-posedness. We pro¬ 
pose finite element methods for the given formulation and prove a priori error 
estimates in Section 3. A numerical example with discretisation errors are 
presented in Section 4 and a short conclusion is drawn in Section 5. 


2 A New Formulation of Poisson equation 


In this section we introduce a new minimisation principle of the Poisson prob¬ 
lem. Let n C d G {2, 3}, be a bounded convex domain with polygonal or 
polyhedral boundary cAI with the outward pointing normal n on dVt. 

We start with the following minimisation problem for the Poisson problem 

Problem 1. Given f G we want to find 


with 


where 


u = argminiL(u) 

v£H^{Q) 


-^(^) = 1 f |Vup dx — ^(u). 
2 Jn 

i{v) = f fvdx. 

Jn 


( 1 ) 

( 2 ) 


We refer to the following references [Tl[2l[6lM IT^[T3l[T^IT6] for different varia¬ 
tional formulations of Poisson equation. 
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Let V = iLo(^i) and Q = and for two vector-valued functions a : 

n —)■ and /3 : O —)■ the Sobolev inner product on the Sobolev space 

iL^(n) {k G M) be defined as 

d 

{o!;P)k,n ■= f3i)k,n, 

i=l 

where (a)* = {/3)i = (3i with a*, /?* G for i = 1, • • • ,d, and the norm 

II ■ ll_f/''(n) is induced from this inner product. We will use the standard notation 
II • lU.n for the norm in the Lf^(fl)-space. We now introduce a functional 
Ja,'r{v,T] /) with 

Ja,-({v, r; /) = ||t||o,q + ||t - aVu||o,Q + j£{v), 

where a > 0 and 7 are two fixed constants, and consider another minimisation 
problem for two variables {v, r) G [IL x Q] 

argmin r;/). (3) 

(D,T)e[y xQ] 


The minimisation problem is equivalent to finding (u, cr) £ [V x Q] such that 
a{{u,a),{v,T)) =-^i{v), (u,T)G[yxQ], (4) 

where the bilinear form a(-, •) is dehned as 

a{{u, a), {v, r)) = (<t, t)o,o + (ct - aVn, r - aVu)o,n- 

Standard arguments can be used to show the continuity of the bilinear form 
a(-, •) on the space F x Q. Now we show that the bilinear form a(-, •) is coercive 
on y X Q. 

Lemma 1. Let a and 7 be two constants with a > 0 . For {u,a) G [IL x Q] 
the bilinear from a(-,-) satisfies 

a{{u,a),{u,(T)) > —(||^x||? f^-F ||cr||2 ) , 

where Ci is the constant in the Poincare inequality 

ll'*^lli,n — C'i||Vu||q^q. 
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Proof. The proof follows from a triangle inequality and Poincare inequality: 


klli,n + ll<’’llo,o ^ + ll<^llo,n 


< 


< 


a 

2^1 pii 

- k 

a 

2Ci + a 
a 

2Ci +a 
a 


— aVullo^Q + IIctIIq^q] + ||cr||Q^f^ 
(lkllo,o + lk-«V^llo,n) 
a{{u,a), {u,a)). 


□ 


Corollary 1. Since the bilinear form a(-,-) is continuous and coercive on 
V X Q, and the linear form i{v) is also continuous on V for f E the 

problem of finding (n, cr) E [P x Q] such that 

a{{u,cr),{v,T)) =(?;, r) E [P x Q], (5) 

has a unique solution from Lax-Milgram lemma. 

Remark 1. In contrast to the standard least-squares method, where we need 
to have I E L‘^{Q), we have here f E Thus the standard least-squares 

method cannot handle the situation if the source function is not , whereas 
the new approach requires exactly the same regularity for f as the standard 
Galerkin approach. 


Let {ue,cre) E P X Q be the solution of the minimisation problem ([3|). We 
now choose a and 7 in such a way that the solution (u, <t) of the minimisation 
problem ([3|). satisfies u = Ug and = Vn. Here the natural norm for an 

element (u, r) E P x Q of the product space P x Q is + Iloilo o- Thus 

([3]) leads to the problem of finding {u, cr) E [P x Q] such that 

2(o-, t)o,o + 2((t - aVu, r - aVu)o,n + J^iv) = 0, {v, r) E [P x Q], 

Letting the test functions r = 0 and u = 0 successively in the above equation 
leads to 

-2{cr - aVu,aVv)o^n-\-'yi{v) = 0, u E P, .g. 

(o-, T)o,n + (cr - oVu, t)o,o =0, r E Q. 

The second equation immediately yields 

(2cr - aVn, T)o,n = 0, r E Q, 
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and hence a = 2 ensures that cr = Vu. Using cr = Vu in the first equation of 
(j 6 |l. we have 

-2a(l — Q:)(Vn, Vu)o,n + ^t{v) = 0. 

We have the standard variational problem for the Poisson equation if 7 = 
2q:(1 — a), and thus setting a = 2, we get 7 = —4. Now we have the following 
problem. 

Problem 2. Given f G H~^{Q), the variational equation for the minimisation 
problem is to find {u, cr) G [U x Q] such that 

a{{u,a),{v,T)) = 2e{v), ('(;,t)g[UxQ], (7) 

where the bilinear form a(-, •) is defined as 

a{{u, cr), {v, r)) = (cr, t)o,o + (cr - 2Vu, r - 2Vu)o,o. 


From the above discussion we have the following proposition. 

Proposition 1. Let u be the solution of Problem\^ and of Problem\^ 

Then we have u = u and a = Vu. 

Remark 2. The idea can be easily generalised to a general differential equa¬ 
tion, which can be put in a minimisation framework. For example, consider 
the solution of the linear elastic problem of finding the displacement field 
u G V = such that 

u = arg min - f e(v) : Ce(v) dx — ^(v), 
vev 2 Jq 

where e(v) = ^ (Vv + [Vv]"^) is the symmetric part of the gradient, C is the 
Hooke’s tensor, and i{-) is a linear form 

£(v) = f i ■ vdx. 

Jn 

By defining a pseudo-stress cr = \/Ce(v), we can put this in the above frame¬ 
work with 

a(u, cr, V, r) = (cr, t)o,o + (cr - 2VCe(u), r - 2VCe(v))o,o. 

We note that since C is a symmetric positive definite tensor, its square is 
well-defined. 
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3 Finite element approximation and a priori error 
estimate 


Let 7ft be a quasi-uniform partition of the domain O in simplices, convex 
quadrilaterals or hexahedra having the mesh-size h. Let T be a reference 
simplex or square or cube, where the reference simplex is dehned as 

d 

T := {x E : Xj > 0, i = 1, - ■ ■ ,d, and Xi < 1}, 

i=l 

and the reference square or cube T := (0, l)*^. 

The hnite element space is defined by affine maps Ft from a reference element 
T to a physical element T £ Th- For k > 0, let Qk{T) be the space of 
polynomials of degree less than or equal to A: in T in the variables xi, • • • , x^ if 
T is the reference simplex, the space of polynomials in T of degree less than or 
equal to k with respect to each variable xi, • • • , x^ if T is the reference square 
or cube. 

Then the finite element space based on the mesh 7ft is dehned as the space of 
continuous functions whose restrictions to an element T are obtained by maps 
of given polynomial functions from the reference element; that is, 

Sh ■= ^Vh E : Vh\T = Lft o E Qk{T), T E Tftj , 

see [^[711101114] , We now dehne 14 := S'ft n We also want to dehne two 

other hnite element spaces Lft and Qft as 

Lft := {vh E L^{n) : v^It G Qk{T), T E Tft} , Qft := [L^f, 

Now a discrete formulation of our problem is to hnd (ttft, cTft) E [14 x Qft] such 
that 

a((uft, cTft), (uft, Tft)) = 2i{vh), (uft, Tft) E [14 x Qft], (8) 

Since [14]“^ C Q, we can also use Vft = [ 14 ]'^ to discretize the gradient of the 
continuous problem. This leads to a problem of hnding (ttft,crft) E [14 x Vft] 
such that 


o((uft,crft), (uft,Tft)) = 2£(uft), (uft,Tft) E [14 X Vft], (9) 

which utilizes equal order interpolation. Since the discrete formulation is con¬ 
forming, the bilinear form a(-, •) and the linear form i(-) are both continuous 
on the corresponding spaces. The coercivity also follows from the continuous 
setting. 
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Theorem 1. Thus the discrete problem of finding {uh,crh) € [Vh x Qh] or 
{uh,cTh) G [Vh X Vh] such that 

a{{uh,(Th),ivh,Th)) = 2£{vh), {vh,Th) ^[VhxQh] or {vh,Th) e [VhxVh] 
has a unique solution, and the solution satisfies 

||n-'W/i||i,o + ||cr-cr/j||o,n < c ( inf ||n - + inf lk - r/i||o,o ) , 

\Vh&Vh Th^'Vh orTh^Clh ) 

where u is the exact solution of the problem © and a = Vu. 

Proof. The proof follows from Galerkin orthogonality and standard arguments. 

□ 

Remark 3. The solution u is assumed to be in Hq{Q) only for the purpose 
of simplicity. In fact, any non-zero Dirichlet condition or mixture of Dirichlet 
and Neumann boundary conditions are all fine as in the case of the standard 
Galerkin finite element method. 


4 Numerical example 


In this section we consider a numerical example to demonstrate the perfor¬ 
mance of this new minimisation scheme. In fact, we show the discretisation 
errors for the solution u in the and R^-norms, and discretisation errors for 
the gradient in the L^-norm. For this example, we consider the domain of the 
square O = [—1,1]^ with the exact solution 

u{x, y) = {{x — y) exp(—5.0(x — 0.5)(a: — 0.5) — 5.0(y — 0.5)(y — 0.5))), 

where the right-hand side function / and the Dirichlet boundary condition 
on dfl is obtained by using this exact solution. The two components of the 
gradient are denoted by and respectively, where their numerical ap¬ 
proximations are denoted by and respectively. We start with the initial 
uniform triangulation of 32 triangles in the first level and then refine uniformly 
in each level. We have tabulated the discretisation errors in Table [TJ We can 
see that the numerical results are as predicted by the theory. Moreover, the 
discretisation errors show the superiority of the scheme as the discretisation 
errors for the gradient of the solution converge quadratically to the exact so¬ 
lution. This is normally not achieved in any mixed finite element methods. 
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Table 1: Discretisation errors for the solution and gradient 


D 

11^.4 - Uh\\i,n 

||w - 



1 

3.41208e-01 


3.71839e-02 


1.72948e-01 


1.72948e-01 


H 

1.70261e-01 

1.00 

1.15857e-02 

1.68 

5.39510e-02 

1.68 

5.39510e-02 

1.68 

H 

8.40661e-02 

1.02 

3.09011e-03 

1.91 

1.45061e-02 

1.89 

1.45061e-02 

1.89 

Q 

4.18485e-02 

1.01 

7.86293e-04 

1.97 

3.74056e-03 

1.96 

3.74056e-03 

1.96 

B 

2.08998e-02 

1.00 

1.97503e-04 

1.99 

9.62306e-04 

1.96 

9.62306e-04 

1.96 

B 

1.04469e-02 

1.00 

4.94393e-05 

2.00 

2.51753e-04 

1.93 

2.51753e-04 

1.93 


5 Conclusion 


We have proposed a new minimisation principle for the Poisson equation based 
on the solution and its gradient. One big advantage of this formulation is that 
a finite element approximation can be performed as in a least-squares finite 
element method without fulfilling the compatibility condition between two 
finite element spaces. However, the finite element approach is much easier than 
in a least-squares approach. An optimal a priori error estimate is proved for 
the proposed formulation. A numerical example is presented to demonstrate 
the optimality of the scheme. 
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